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Technical Objectives 

The Hackett method-a Wall Pressure Signature method has been selected to be adapted 
for the 12-ft Wind Tunnel WIAC system in the present phase. This method uses limited 
measurements of the static pressure at the wall, in conjunction with the solid wall 
boundary condition, to determine the strength and distribution of singularities representing 
the test article. The singularities are used in turn for estimating wall interference at the 
model location. Hackett’s method will have to be formulated for application to the unique 
geometry of the 12-ft Tunnel. The development and implementation of a working 
prototype will be completed, delivered and documented with a software manual in this 
phase. 

The WIAC code will be validated by conducting numerically simulated experiments rather 
actual wind tunnel experiments. It is an effective, but efficient way to validate the 
implemented code. The simulations will be used to generate both free-air and confined 
wind-tunnel flow fields for each of the test articles over a range of test configurations. 
Specifically, the pressure signature at the test section wall will be computed for the 
confined case to provide the simulated "measured" data. These data will serve as the 
input for the WIAC method-Hackett’s method in this phase. The performance of the 
WIAC method then may be evaluated by comparing the corrected parameters with those 
for the free-air simulation. Each set of wind tunnel/test article numerical simulations 
provides data to validate the WIAC method. 
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Status of Progress 

1. System Validation Design 

A numerical wind tunnel test simulation is planned to be utilized to validate the WIAC 
methods developed in the project, specifically to the wall pressure signature method-- 
Hackett’s method in the current phase. A low-order potential-flow panel method has been 
selected to simulate the wind tunnel and test article geometries for its flexibility. 

A copy of the low-order panel code, called PMARC (Panel Method Ames Research 
Center), was obtained from Dr. Frank Steinle, technical monitor. The code of PMARC 
(Ref. 1) is a well-documented code with open architecture which has several advanced 
features including internal flow modeling for wind tunnel application. The code has 
another feature allowing the adjustment of the size of panel for suiting the computer 
hardware and the size of the problem. PMARC can be run on computers ranging from 
a Macintosh II workstation to a Cray Y-MP. 

Currently, the copy of PMARC (version 12) has been compiled on a IBM R6000 and 
checked out on two sample cases -- wing and wing/body which have been provided by 
NASA/ARC. The PMAPP, a post processor for PMARC, will be received from 
NASA/ARC shortly. 

2. Hackett’s Method Implementation 

The implementation of a Wall Pressure Signature method which is based upon the wall 
interference method developed by Hackett (Ref. 2) has been initiated. In the present 
reported period, the blockage correction is developed and implemented. The lifting 
correction will be carried out in the next period. A working computer code has been 
implemented and verified by an analytic solution. 

The analytic flow field solution of a rectangular wind tunnel with a finite airfoil wing and 
its wake was derived as a bench mark solution to verify the solution calculated by 
Hackett’s Method. The blockage correction of the Wall Pressure Signature method 
(Hackett’s Method) for this wing/wake case is obtained and verified by the analytical 
solution. 

The detailed development and results are given in Attachment I. 


Future Plan 

1 . Panel Code Application 
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The application of panel code PMARC will be initiated for Ames 12-ft Wind Tunnel. The 
wall pressure in using Hackett’s method will be computed for the cases of wing and 
wing/body configuration. The post processor (PMAPP) for PMARC will be installed on 
the UTSI’s computer to support PMARC’s application. The software system validation 
will be started for Hackett’s method developed in Task 2 with the configuration of 12-ft 
Tunnel. 

2. Lifting Correction of Hackett’s Method 

The development and implementation of the lifting Correction of Hackett’s Method will be 
initiated. The verification will be based on an independent solution such as a 
corresponding solution of lifting case derived in this reported period for the blockage case. 
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BLOCKAGE CORRECTION IN 3-DIMENSIONAL WIND TUNNEL 
TESTING BASED ON THE WALL SIGNATURE METHOD 
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BLOCKAGE CORRECTION IN 3-DIMENSIONAL WIND TUNNEL 
TESTING BASED ON THE WALL SIGNATURE METHOD 


N. Ulbrich and C. F. Lo 
The University of Tennessee Space Institute 
Tullahoma, Tennessee 37388 

Summary 

A blockage correction technique in 3-dimensional wind tunnel flow fields is developed 
and verified. Blockage corrections are obtained by Hackett’s method with pressure mea- 
surements on the wind tunnel wall. 

An analytic wind tunnel flow field solution of a rectangular wing with parabolic airfoil 
section and a trailing edge wake is presented. This analytic solution is used to simulate 
pressure measurements on selected wall locations and to calculate the exact solution of the 
blockage problem. 

A pressure coefficent correction due to model blockage is obtained by Hackett’s method 
in using simulated wall pressure measurements. Pressure coefficient corrections show good 
agreement with the exact solution of the blockage correction problem. This verifies the 
implementation of Hackett’s method for the blockage correction. 

1. Introduction 


The Wall Pressure Signature Technique (Hackett’s method) 1, 2 ’ 3,4, 5 is considered to 
calculate blockage correction due to wall interference in the NASA/ARC 12ft Pressure 
Wind Tunnel. The first part of this report presents the development and implementation 
of an analytical correction technique which requires the simulation of wall pressure mea- 
surements for Hackett’s method. A wind tunnel with rectangular cross section is selected. 
A rectangular wing with parabolic airfoil section is chosen as a test article and a solution 
of this flow field is obtained by applying thin wing theory. The influence of the wake is 
modeled by using a line source and a line sink. The analytic solution of the wind tunnel 
flow field of a wing and a wake is obtained by combining free-air solutions with the method 
of images. 

Two important sets of data are derived from the analytic solution: wall pressure mea- 
surements are simulated which will be used as required by Hackett’s method to calculate 
blockage corrections; the exact solution of the blockage problem is calculated to verify 
Hackett’s method. Basic steps of the development and verification of a blockage correction 
technique based on wall pressure measurements are shown in Fig. 1. 

A detailed description of Hackett’s method is presented in the second part of this 
report. It is shown how an equivalent flow field representation is defined based on wall 
pressure measurements. This equivalent flow field model is then used in combination 
with the method of images to calculate a blockage correction. This blockage correction is 
compared to the exact solution of the blockage problem to verify Hackett’s method. 
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2. Analytic Wind Tunnel Flow Field 

2.1 General Remarks 


The simulation of a three dimensional wind tunnel flow field requires the mathemati- 
cal description of a simplified test article and a wake. A rectangular wing with parabolic 
airfoil section is selected as a test article. In general the analytic solution of a wing in in- 
compressible three-dimensional flow cannot be obtained easily. Integration of complicated 
algebraic expressions is necessary. 

The solution of the thickness problem of a thin wing in incompressible free-air flow is 
known in the form of a double integral. Basic equations will be presented and the geometry 
of the selected rectangular tunnel and wing will be specified. The free air solution is 
obtained by direct integration. A wake in free-air flow can be simulated by using a finite 
span line source and line sink. An analytic flow field solution will be derived by direct 
integration as well. 

The wind tunnel flow field will be calculated by superimposing the analytic solutions 
of wing and wake and applying the method of images to the total flow field to satisfy the 
wall boundary conditions. 

Finally formulae of the analytic solution of the blockage problem are obtained. 

2.2 Free-Air Flow Field of the Wing 


The solution of the perturbation potential of three-dimensional planar wings in in- 
compressible free-air flow is given by Ashley and Landahfi as follows 




J_ff dg(ri,yi) 

2 7T JJ dxi 

A 


dx i dy i 

\J{x~— ij ) 2 + (y — yi ) 2 + z2 


( 1 ) 


where ‘A’ is the projected wing area, ‘r’ is the thickness and dg{x \ , y\ )/ dx\ is the thickness 
slope. 

Figure 2 shows a rectangular wing with parabolic airfoil section, wing span ‘s’ and 
chord length l c\ The thickness function y(ii,yi) of the parabolic airfoil section has the 
following form 


y(* i>yi) 

The thickness slope is then given as 


2 • ij 
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Combining Eqs.(l), (2b) and knowing that 0 < xj < c and 


2x\ 
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The dimensionless axial perturbation velocity is given as 


uj(x,y,z) 


dfi 

dx 


r f f [1 - 2xi/c] ■ (x - Xi) 
n j _j^[(x-xi ) 2 + (y-yi) 2 +z 2 ] 3/2 


( 4 ) 


or in terms of the corresponding pressure coefficient 

c Pl (x,y,z ) = — 2 • uj(x,y,z) (5) 

The integration of Eq.(4) is complicated. Different types of integrals have to be 
considered depending on the choice of field point. 

In general the solution consists of six cases to cover the entire flow field. The classi- 
fication of these cases and the solution of the incompressible axial perturbation velocity 
are given in Appendix A. Singularities exist at the leading and trailing edge of the wing. 
Otherwise a well defined solution can be found in the flow field. 

A verification of parts of the solution is possible if the limit as the wing span goes to 
infinity is taken. The axial perturbation velocity on the model surface is obtained based 
on Eq. (A. 3) and limits given in Appendix B. We get 

u(x,y,0) = |-r.(l + \ • [1 - 2i/c] - In jj^L) (6) 

Equation (6) agrees with the solution of a two-dimensional parabolic airfoil section given 
by Schlichting and Truckenbrodt 7 . 

Pressure coefficients in free- air on selected locations are calculated based on Eq. (5) 
and the analytic solution given in the Appendix A . The following wing geometry is chosen: 
chord length c = 1.0, wing span s = 4.0, thickness r = 0.1. Pressure coefficients are 
calculated on the planes z = 0.0, 0.2. Results are presented in Figs. (3a) and (3b). 

Singularities on the leading and trailing edge of the wing can be noticed in Fig. (3a). It 
is obvious that the pressure coefficient distribution is symmetric with respect to the planes 
y = 0.0 and x = c/2. This symmetry property can be used to reduce calculation time if 
pressure coefficients in a plane with constant ^-coordinate are determined. The absolute 
value of the pressure coefficient decreases with inceasing z-coordinate as the perturbations 
of the free-air flow field decay with increasing distance from the wing plane. 

2.3 Free- Air Flow Field of the Trailing Edge Wake 

The flow field of the trailing edge wake of a rectangular wing is simulated using a line 
source located on the trailing edge and a line sink of equal strength at a large distance ‘ c 9 ’ 
downstream of the trailing edge. The geometry of the wake model is given in Fig. (4). 

The general solution of the velocity potential of a single line source of span 4 b and 
constant source strength 4 g’ in free-air flow is given by Katz and Plotkin 8 as 
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<p(x,y,z) = - 


*2 

I \ Ax' 


di 


4?r J yfix - *i ) 2 + (y - yi ) 2 + (* - 2 i ) 2 

fi 


(7a) 


We now consider a line source and a line sink located in the x-y plane as shown in Fig. 4. 
In this case we get for the potential the following equation 
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The axial perturbation velocity U 2 (x,y,z) can now be obtained by diffentiating potential 
<P2 with respect to x. We get then 


u 2 (x,y,z) 


&£2 

dx 


(8a) 


or 


V /2 
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y+b w /2 


[(i - x a ) 2 + z 2 + i 2 f /2 


y—b w /2 


( 86 ) 




di 




[(x-x a -c a ) 2 + 2 2 + £ 2 ] 3/2 


The complete analytic solution of the free-air solution of a wake is given in Appendix C. 
The source/sink strength l q ’ is equal to the wake thickness *d w ' if is dimensionless. 


2.4 Wall Pressure Measurements and Blockag e 


The solution of the axial perturbation velocity of a rectangular wing and a trailing edge 
wake in a rectangular wind tunnel of width W and height H (see Fig. 2) can be derived 
if the method of images is combined with the free-air solution. The tunnel solution is 
obtained by superimposing free-air solutions as follows: 


OO OO 

ut(x, y, z) = E E uj(x,y 4- mW,z + nH) 

m= — oo n = — oo 

OO oo 

+ u 2 (x,y + mW,z + nH) 

m = — oo n= — oo 


(9a) 
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where the velocity ‘ui’ is given by Eqs.(A.l) ... (A.5) in Appendix A and velocity ‘u 2 ’ 
is given by Eqs.(C.l) and(C.2) in Appendix C. It is only necessary to substitute ‘y’ by 
‘y + mW' and V by l z + nH' in these equations. 

For practical calculations Eq. (9a) is approximated as follows 

k k 

ut(x, y, z) « ]|P ui(x ,y + mW,z + nH) 

m= — k n=—k 

k k 

+ E E u 2 (x,y + mW,z + nH) 

m= — k n= — k 

where ‘fc’ has to be a large number to guarantee convergence of the solution. The corre- 
sponding pressure coefficient in the flow field is given as 

c PT (x,y,z) = -2 -u T (x,y,z) (9c) 

A computer cannot represent numbers exactly and errors are introduced due to finite 
digit arithmetic. The wind tunnel flow field of the wing is considered alone to study this 
problem. The wing geometry of the free-air study, a tunnel width W = 6.0 and a tunnel 
height H = 4.0 are selected. A single point in the flow field is chosen ( x = 0.5, y = 0.0, 
z = 0.2) and the pressure coefficient is plotted as a function of the chosen summation limit 
‘A:’ (see Fig. (5)). The single precision solution becomes unstable at k m 8. The double 
precision solution converges uniformly. 

Figure 6 shows a specific wind tunnel, wake and wing geometry. For this case pressure 
coefficients are calculated in the plane z = 1.75 (upper wind tunnel wall). Results are 
presented in Fig. 7. It is now possible to simulate input data for the application of Hackett’s 
method if a single row on the tunnel wall is selected (e.g. y = 0.0 and z = 1.75). Wall 
pressure coefficients are calculated on 31 different locations on this row. Calculated results 
are depicted as the solid line in Fig. 8. 

The verification of blockage corrections obtained by Hackett’s method requires a com- 
parison with the exact solution of the blockage problem. The exact solution of the blockage 
problem of the given superimposed flow field can be calculated based on the free-air and 
wind tunnel solutions. This correction on the plane z = 0 is defined as follows 

c P( .(x,y,0) = c pr (x,y, 0) - c p (x,y,0) (10a) 

where c PT is given by Eq. (9c). The free-air pressure coefficient is calculated based on 
Eqs. (4) and (8b) as follows 

c p = - 2 • [ uj + u 2 ] (106) 

The calcuation of the pressure coefficient correction has to be done in double precision 
arithmetic to achieve sufficient accuracy. 

The blockage correction c pi is now calculated for a single row on the surface of the 
model (y = 0). These results are depicted as the solid line in Fig. (9). This analytic 
solution of the blockage problem will be used to verify Hackett’s method. 


(96) 
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3. Hackett’s Method 

3.1 General Remarks 


Wall pressure signature methods like Hackett’s method were successfully applied in 
the past to calculate blockage corrections in 3-dimensional wind tunnel testing. Hackett’s 
method consists of two major parts: the signature analysis of the wall pressure mea- 
surements to model test articles and the calculation of blockage corrections based on the 
simplified test article representation obtained from the signature analysis. 

At first a description of the signature analysis will be presented. The signature analysis 
is applied to a simulated wall signature which is calculated for a given wing, wake and wind 
tunnel geometry as shown in Fig. 6. Secondly the simplified model/wake representation 
is derived. Finally, blockage corrections are calculated based on the simplified model 
representation and the method of images which are compared to the exact solution of the 
blockage problem. 

3.2 Signature Analysis 

These simulated measurements of a selected model/wake/wind tunnel geometry are 
obtained from the analytic solution described in Section 2.4 . Simulated wall pressure 
coefficients on a selected wind tunnel wall (i.e. —2.5 < y < 2.5 and z = 1.75) are given in 
Fig. 7. The signature analysis suggested by Hacketf is applied to simulated wall pressure 
measurements on a single row (in our case y = 0.0 and z = 1.75; see Fig. 6). 

The wall pressure measurement on the selected row (solid line in Fig. 8) is split into 
a symmetric and antisymmetric part by iteration. This iteration requires that the peak 
location of the symmetric signal is identical with the location of the inflection point of the 
antisymmetric signal. 

The symmetric part of the signal is related to a simplified representation of the solid 
body blockage in the form of a line source/sink pair. The source/sink spacing and location 
and the source sink strength can be determined from the width at half height of the 
symmetric signal, the location of the velocity peak and the maximum value of the velocity. 

The antisymmetric part is similarly related to a simplified representation of the wake. 
Characteristic parameters (like line source/sink location and source/sink strength) are de- 
rived from the location of the inflection point of the antisymmetric signal and the asymp- 
totic value of the velocity far downstream of the model. More details on the signature 
analysis are given by Hackett 2 . A concise description of the method is also given by 
Allmarai i 5 . 

3.3 Blockage Correction 

Blockage corrections sire now calculated based on the simplified test article and wake 
representation which is given in terms of 4 line sources /sinks. The equations of the blockage 
correction of the simplified test article and wake representation can be derived easily by 
modifying equations presented in Appendix C and using the method of images. The result 
of this calculation is given by the dashed line in Fig. 9 and compares well to the known 
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analytic solution of the blockage problem. Pressure coefficient corrections given in Fig. 9 
have to be combined with wind tunnel measurements on the same location in the x-y-z 
coordinate system to obtain the free-air measurements. 

4. Remarks 

Blockage corrections in a 3-dimensional wind tunnel flow field axe calculated applying 
Hackett’s method to simulated wall pressure measurements. These corrections show excel- 
lent agreement with the known analytic solution of the blockage problem. The influence 
of a fuselage on the applicability of Hackett’s method still has to be included in these pre- 
liminary results. An improvement of the signature analysis is also considered to simplify 
the calculation of the equivalent test article and wake representation. 
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Fig. 1 Development and Verification of a 3-Diinensional 
Blockage Correction Technique 
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APPENDIX A 


ANALYTIC FLOW FIELD SOLUTION (WING) 


Case No. 

Condition 

Solution 

1 

x arbitrary 

\lj\ 7^ 5 / 2 
2^0 

Eq. (A.l) 

2 

x arbitrary 

Itrl = */2 
2 0 

Eq. (A. 2) 

3 

i/O and r/c 

M 7 ^ s/2 
2 = 0 

Eq. (A. 3) 

4 

i/O and x ^ c 

|y| = s/2 
2 = 0 

Eq. (A.4) 

5 

x = 0 or x = c 

|y| > «/2 
2 = 0 

Eq. (A.5) 

6 

I 

i = 0 or i = c 

Ivl < s/2 
2 = 0 

singular 
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APPENDIX B 


CALCULATED LIMITS 
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= 1 


lim 

§—►00 


y/x 2 + (y ± s/2) 2 - |y ± s/2| 

[ ^/(x - c) 2 + (y ± s/2) 2 - |y ± s/2|J 


(x - c) 


lim 


(y + s/2) ■ In 


x + ^x 2 + (y + -s/2) 2 


_(x - c) + y(x - c) 2 + (y + s/2) 2 _ 


lim 

§—►00 


(» - s/2) • In 


i + + (y - s/2) 1 


.(* - c) + \/(x - c) ! +(v - s/2) 3 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B. 5 ) 
(B.6) 
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APPENDIX C 


ANALYTIC FLOW FIELD SOLUTION (WAKE) 


A, = 


Case No. 

Condition 

Solution 

1 

z ^ 0 or 
2 = 0 and 

i ^ i, and x ^ x t + c. 

Eq. (C.l) 

2 

2 = 0 and 

x = x s or x = x 3 + c 9 and 

|y| > b w /2 

Eq. (C.2) 

3 

z ~ 0 and 

x — x g or x — x g “f* c g and 

|y| < b w /2 

singular 

= 

— A'l ] ! 

y + b w /2 

(C.l) 

-x a ) 2 + 2 2 

y - 

yj{x - x,) 2 + (y + b w /2) 2 + z 2 
b w /2 ] 


\J{x - x,) 2 + (y - b w /2) 2 + z 2 


A<i = 


y + K! 2 


(x-X t -C,) 7 +Z* + (y + bj 2f + Z 2 

y - b w /2 

\[{x - ~x s - c,) 2 + (y - &u,/2) 2 + 2 2 J 


= 4r£ 


y + W 2 


y — ^u»/2 


\/c, 2 + (y + &„/2) 2 


(C.2) 
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